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Abstract 

We consider a hierarchical nested microstructure, which also contains a point of sin¬ 
gularity (disclination) at the origin, observed in lead orthovanadate. We show how to 
exactly compute the energy cost and associated displacement field within linearized elas¬ 
ticity by enforcing geometric compatibility of strains across interfaces of the three-phase 
mixture of distortions (variants) in the microstructure. We prove that the mechanical 
deformation is purely elastic and discuss the behavior of the system close to the origin. 

Keywords. A. microstructures, phase transformation. B. strain compatibility. C. asymp¬ 
totic analysis, variational calculus. 
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1 Introduction 

Solid-to-solid phase transformations are often accompanied by the formation of unusual and 
intriguing mixtures of phases at the mesoscale spanning nanometers to microns in length 
scales [7]. In the case of metallic alloys, below the transition temperature it is common to 
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observe the coexistence of fine layers of martensitic or product twins (or ’’variants”) with the 
parent austenite phase of higher symmetry. Martensitic transformations are displacive, are 
driven by shear strain and/or shuffles (intracell atomic displacements), and are invariably 
first-order in nature that leads to hysteresis and metastability [3]. The rich microstructure 
seen in high resolution electron microscopy (HREM) is often a manifestation of this metasta¬ 
bility [iniis]. Various approaches have been utilized over the last several decades to model the 
emergence of microstructure in martensites. They are largely variational using a free energy 
potential, and either use finite deformation, sharp interfaces and iteratively minimize a free 
energy, or start with random initial conditions and evolve a free energy potential according 
to some dynamics m- The traditional phase field approach is often used in the limit of small 
strains, and methods based on Ginzburg-Landau theory can be applied for small strains as 
well as finite deformation. 

The study of mixtures in the framework of a variational setting traces back to the work 
of Ball and James [1]. Their technique consists of matching different crystal phases, possibly 
at different scales, through geometric compatibility. In the literature, this idea has been 
widely applied in the analysis of periodic mixtures, the situation in which the physical and 
geometrical properties of these mixtures are repeated periodically in an elastic body. 

Even though the case of periodic microstructure is of central importance in the modeling 
of composite materials, and it has become a classic subject of study, the situation related to 
more general microstructure, and specifically the case of self-similar fine hierarchies, remains 
to be fully explored. In fact, there are a number of fascinating open problems and questions, 
among which is the study of a large family of fine hierarchical structures observed in artificial 
polymers and biological materials (e.g., bones and leaves). Our emphasis will be on the family 
of heterogeneities in which there is an interaction of topological singularities that leads to 
fascinating (non-periodic) microstructure. 

In this paper we focus on the class of hexagonal-to-orthorhombic transformations where 
three equivalent stretching directions of the parent austenite phase give rise to three orien¬ 
tations of the product martensite. Examples of compounds undergoing this transformation 
include the mineral Mg-cordierite, Mg 2 Al 4 Si 50 i 8 , and Mg-Cd alloys [S]. Closely related 
materials include those undergoing a hexagonal-to-monoclinic transformation, such as lead 
orthovanadate, Pb 3 (V 04 ) 2 , and samarium sesquioxide, Sm 203 [10]. As the variants need 
to rotate to match at the domain walls, and the domain walls connecting the variants may 
intersect, these materials provide us with an excellent opportunity to study disclinations in 
crystals. Disclinations are formed when the nodes generated by the intersection of the domain 
walls do not close to an angle of 2tt. 

We will study the hexagonal-to-orthorhombic transformation in two dimensions (2D) for 
which the corresponding transformation is triangle-to-centered-rectangle. As the compounds 
consist of stacking or layering of tetrahedral units, the microstructure is essentially homo¬ 
geneous perpendicular to the plane of the paper and therefore 2D is justified. One of the 
most intriguing microstructures observed in this transformation is the self-similarly nested 
tripole-star pattern (see Fig. [^(a)). These transformations have recently been the object 
of extensive numerical study. In Ref. [12] the modeling is based on the minimization of a 
non-convex Ginzburg-Landau potential, both in the scenario of finite and infinitesimal elas¬ 
ticity. It can be observed that in this microstructure the deformation gradient is (nearly) 
piecewise constant, with three possible strain values that correspond to the three wells of 
the free energy density. As a first approximation, the sets where the deformation gradient 
is constant are particular kyte-shaped polygons. These polygons are all identical up to a 
rotation and a rescaling, close to the center of the star, and their measure tends to zero. 

A thorough theory that treats microstructure and other phenomena, such as disclinations, 
dislocations, cavitations or cracks simultaneously, is to the best of our knowledge, still lack- 
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ing. Such a theory would provide an important tool to study various classes of phenomena in 
materials science. Focusing on the case of linearized elasticity, the main result of this paper 
(Theorem is the exact computation of the zero-energy microstructure and of the displace¬ 
ment field that realizes it. In our construction we adapt the techniques of Ball and James 
[U (see also m) to the case of a non-periodic microstructure with a laborious construction 
from scratch. We note that the situation presented is a generalization of the case of simple 
laminates (matching of martensitic variants at the microscopic scale) or possibly, to the case 
of laminates-within-laminates. Indeed, the nesting of the microstructure requires a slightly 
more delicate construction consisting of simultaneous matching of phases across interfaces 
of polygons of different size at an infinite number of scales. Furthermore, matching of ge¬ 
ometrically compatible variants is typically achieved by piecewise affine displacement fields 
(laminates). Continuity of the displacement field is a key property of elastic models because 
it rules out the occurrence of irreversible phenomena such as cavitations and cracks. In the 
present scenario, whether a continuous displacement can realize the microstructure depicted 
in Figj^(a) according to our model, is not a priori clear. The reason is that the center of 
the star may act as a topological singularity for the microstructure. A by-product of our 
analysis in Theorem is that the displacement u is indeed continuous and, consequently, it 
is a purely elastic displacement. 

Whether adopting a geometrically linearized model is physically sound for studying this 
type of microstructure is discussed in Ref. [12] and in the final section of this paper. Al¬ 
though models in linearized elasticity have some intrinsic limitations that may rule out the 
understanding of irreversible or higher order phenomena, we remark here that the intrinsic 
simplicity of this model may shed some light on certain features that survive upon lineariza¬ 
tion, such as the geometry of the microstructure. The analysis of the full non-linear model 
is left out from this article and is the object of ongoing research. 

The paper is organized as follows. In the next section we introduce the basic notation 
used through the paper and the basic concepts of finite elasticity. Moreover, we present 
the mechanical model used in Ref. [12] and its geometrical linearization obtained in the 
asymptotic expansion of small displacements. In Section [^ we present the full analysis of 
the geometry of the problem and the construction of a possible displacement field u which 
reproduces a microstructure similar to the one observed in Fig. [^(a) (Theorem]^. Finally, 
Section]^ is devoted to the physical interpretation of the results contained in Theoremand 
to the presentation of some open problems. 

2 Background 

2.1 Notation 

We gather here the main symbols and the notation used throughout the paper. Our main 
references are [3] and jl). Let N and K denote the set of natural and real numbers re¬ 
spectively. For any integer n, K" is the space of n-dimensional vectors with canonical 
basis {ii};* = 1,.. .n, with origin O = (0,... ,0) and the space of square real ma¬ 

trices. The determinant, the trace and the transpose of the matrix F in M"^” are de¬ 
noted by detF, trF, F^ respectively. We endow with the usual inner product 

F : M := tr (FM^) = corresponding norm |M| := (M : M)^/^. Here 

Mij,Fij are the cartesian components of M and F. The identity in is denoted with 

I with components Sij. We have the orthogonal decomposition of a matrix F G 
M = Msym -I- Msfceuj where M^ym '■= (M -t- M^)/2 and Mskew '■= (M - M^)/2. By further 
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Figure 1: (a): Experimental observation of the tripole-star microstructure in lead orthovana¬ 
date [ini. The two stars are modeled as rank-one three-phase martensitic mixtures, (b): 
Numerical solution of the mechanical equilibrium equations corresponding to the non-linear 
elastic model defined in Eqs. (2.6) and (2.9). (c): Comparison between the strain fields 
obtained numerically in the approximation of linear elasticity (dotted lines) and with the full 
nonlinear model (solid lines) (see also [H]). The graphic shows the strain profiles along a 
line which crosses the tripole-star microstructure (indicated with a dotted line in (b)). In 
both cases the choice of the parameters of the model yields e « 0.156. 
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( 2 . 1 ) 


decomposition of into its deviatoric and spherical part we have 

M = 'M.dev + Mgp/i + Msfceuj 

where := and M.dev := - M^ph- 

We denote with Bii{0, 0) the disk in with center in the origin and with radius equal 
to R: 

Br{0, 0) := {{x, -.x^ +y'^ < R^}. (2.2) 

Let fl be an open and bounded subset of K". We denote with C(n, K") the space of continuous 
vector fields and with (7(11, M") those maps in (7(H,IR") which are continuous on the closure 
of n. Letting p G [1, +oo), we introduce L^(n), the space of measurable functions u : H i—>■ K 
such that \uY‘dx < +oo. Analogously, LP(r2,IR") and L^(f2, M"^”), respectively the spaces 
of vectors or matrices with components in LP(H). Then, K") is the spaces of vector¬ 

valued L^’-functions whose gradient has LP-integrable components. The space W^’°°(f2, K”) is 
that of vectors with essentially bounded components whose gradient has essentially bounded 
components. For regular enough we identify VF^’°°(f2, K”) with the space of Lipschitz 
functions over SI. Other spaces of functions may be defined when encountered throughout 
the paper. 


2.2 Finite elasticity 


According to Ref. |3] Chapter 2 we introduce the deformation gradient F whose components 
are defined as 

dXr 


F- ■ — 

“ dx,' 


where Xi is the Ah component of the position vector of a mass element in the current 
configuration and Xj is the jth component of its position vector in the initial reference 
configuration. In the case n = 2 we let xi = x, X2 = y- The deformation gradient can be 
written in terms of the displacement gradient as 


F ■ — 


^{{.Xi Xi) -\- Xij d{ui Xi) 


dxj 


dxj 


duj 

dxi 


(2.3) 


where Ui are the components of the displacement vector u. The Lagrangian strain tensor is 
defined as E = 1(F^F — I). In components, the Lagrangian strain tensor reads 


Ifdui duj duk duk \ 

= d ( AW + AW + -’ 


2 V dxj 


dxi 


dxi dxiJ 


(2.4) 


Finally, we present the symmetry-adapted Lagrangian strains expressed, as in m Section 
6.1], as functions of the displacement gradient: 
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^xy + Eya ; _ l/pp dUy 

-2- - 2 + ^21^22) - 2 


dx 


dux dux 
dx dy 


duy duy 
dx dy 


(2.5) 
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2.3 The Ginzburg-Landau model for the triangle-to-centered- 
rectangle transformation 

In this Section we briefly review the Ginzburg-Landau model for the 2D triangle-to-centered- 
rectangle (TR) transformation introduced in Refs. [I2l[6l|9j. The TR transformation is the 
2D version of the 3D hexagonal-to-orthorhombic transformation that occurs in materials such 
as the MgCd ordered alloy. A peculiar aspect of these materials is that they are rare examples 
of disclinations in crystals. In the 2D setting, disclinations are point defects that result when 
multiple twin boundaries intersect. The matching condition between two martensite variants 
requires the crystal lattice to rotate, and the total rotation angle along a closed circuit that 
contains the intersection of the twin boundaries (disclination) is, in general, nonzero. This 
requires a stretch of the lattice, additional to the transformation strain, in order to maintain 
its coherency. The disclination is characterized by the disclination angle, which is the total 
rotation of the crystal lattice along the closed circuit that contains the defect. 

A special case of disclination and the most interesting pattern in the microstructure of the 
MgCd alloy is the self-similarly nested tripole-star disclination, also observed in lead ortho¬ 
vanadate, Pb 3 (V 04 ) 2 , undergoing a trigonal-to-monoclinic transformation [TU]. This pattern 
was analyzed in detail by Kitano and Kifune [S]. In this microstructure the intersection of 
twin boundaries do not require a stretch of the lattice for their matching, and therefore they 
are not disclinations. However, the star pattern is itself a disclination, as the total rotation 
angle along a closed circuit containing the center of the star is nonzero. 

In order to study this microstructure, in this paper we adopt the Ginzburg-Landau ap¬ 
proach appropriate for the TR transformation. The corresponding Landau free energy den¬ 
sity, is a function of the symmetry adapted components of the Lagrangian strain tensor, 
and has the functional form 



Figure 2: Level curves of '0^ for ei = 0. Here A = 1, B = —30, C = 200, T = O.STc, Tc = 1 
(e« 0.156). 

62 , 63 ) = -^63 -I- —[T — Tc){e 2 + 63 ) + -^(62 — 36263 ) -I- -(62 + 63 )^ +M. (2.6) 

The phase transition is first order for B ^ 0 and the transition temperature is Tg = Tq + 
2B‘^19AC. For T < Tq the free energy density has three absolute minima, corresponding to 
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the three martensitic variants, with strain values, 


and 


ei =0, 62 = e, 

63 = 0 

Cl =0, 62 = 

63 = ■\/3e2 

ei =0, 62 = 

63 = —■\/3e2, 

-B + - Tc) 

^ 2C 


(2.7) 


For < 0, e is a positive (real) number. The parameter j\4 in (2.6 1 is chosen so that 
min'ipL = 0. Thus, 

M = -['^{T-Tc)e^ + ^e^ + ^e* 

Moreover, variations in the order parameter strain fields are incorporated through a gradient 
term. Thus, the total free energy density considered is 


ipGL = ipL + y (|Ve 2 (a;,y)p + \\/e3{x,y)\‘^) . 


( 2 . 8 ) 


The approach adopted in Ref. [T^ consists of minimizing the total energy of the system. 


= 


/a 


K 

'<pL{ei{x,y),e2{x,y),e3{x,y)) + y (|Ve2(x, y)|^ 


|Ve 3 (a;,y)n dxdy, (2.9) 


with the constraint of compatibility of the strain fields. This can be done using a relaxational 


algorithm or solving the Euler-Lagrange equations associated with (2.91. 

In Figj^(b) we show a snapshot of a microstructure with several tripole-star patterns 
obtained as a solution of the mechanical equilibrium equations using an iterative spectral 
method [12j . The microstructure is obtained with e « 0.156 and only the strain component 63 
is shown, with light (dark) greyscale corresponding to positive (negative) values of 63 . Notice 
that 63 is (approximately) constant over each subregion of the same color. The same behavior 
is also observed for 62 - On the contrary, the strain component ei is (approximately) zero 
within the subregions, whereas nonzero at the boundaries of the subregions (twin boundaries). 
While the effect of the regularizing term ^(|Ve 2 p-l-|Ve 3 p) is to smooth the twin boundaries 
separating the different subregions, the values of the symmetric Lagrangian strain components 
within each subregion are determined by ipL. 

The full non-linear model (Eq. (2.9)) is appropriate for the analysis of disclinations in 


crystals, as it can distinguish the different types of interfaces present in the microstructure 
which are the origin of these defects. Although a geometric linearization of this model 
does not capture the existence of the disclination, the linear model is still able to reproduce 
the tripole-star pattern |12j . In the next Subsection we recall the approximations made in 
geometric linear elasticity and apply them to the Ginzburg-Landau model introduced in this 
Section. 
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2.4 Geometric linearization of the Ginzburg-Landau model 


In the hypothesis of small displacements (and small displacement gradients) one can actually 
neglect the higher order terms in (2.5) yielding 


ei = 

^xx ^yy ^ 

1 

'dUj; 


dUy ■ 

2 

" 2 

dx 


dy \ 

62 = 

^xx ^yy ^ 

1 

-dux 


dUy ■ 

2 

" 2 

dx 


dy J 

63 = 

^xy ^yx 

^ 1 

s 


dUy ■ 

2 

” 2 

L dy 


dx . 


:= 

■= £2{x,y), 


■■= £3{x,y). 


( 2 . 10 ) 


Upon introduction of E = E(Vu) := Vusym the (linearized) mechanical strain tensor, we 
have that £2 and £3 are the components of the deviatoric part of E and £1 is the component 
of the spherical part of E: 


E = 


Vu^) 


£2 £3 

£3 -£2 


+ £il. 


( 2 . 11 ) 


In the linearized elasticity scenario, the new mechanical model is simply obtained by plugging 
£ 1 , £2 and £3 instead of ei, 62 and 63 respectively into the Ginzburg-Landau free energy density 
defined in Eqs. (2.6) and (2.8). The Landau free energy density, which we still denote by 
^z,, reads 

V'l(£ii£ 2,£3) = “ ^c )(£2 + £ 3 ) + ^(£2 “ 3£2£3) -(£3 + £ 3 )^ + Af. (2.12) 


Minimization of (2.12) yields, in terms of strain matrices written as in (2.11), the three 
minimum points 



Notice that even in the linearized model tpL has three distinct minimum points. Importantly: 

|Ei|2 = |E2|2 = |E3p = 2e2. (2.14) 


3 The tripole-star pattern 


Minimization of the total energy 



K 

'^L{£i{x,y),e2{x,y),e3{x,y)) + (|V£2(x, j/)|" 


|V£3(a:,?/)p) dxdy 


(3.1) 


yields numerical solutions for the strain fields which are, at least from a qualitative point 
of view, analogous to those obtained for the full nonlinear model m- Briefly, the term ipL 
enforces the symmetrized gradient of u to be very close to the either Ei,E 2 or E 3 while 
the presence of ^(|V£ 2 p -I- |V£ 3 p) prevents the formation of steep interfaces between the 
three phases. Additional models are available in the literature which penalize the length 
of the interfaces and at the same time allow jump discontinuities of the variants. We refer 
specifically to a forthcoming paper of Dr. A. Ruland for an analysis of a model with inter¬ 
facial energy terms for the tripole-star considered in this paper. In Figj^(c) we compare 
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the strain fields obtained with the linearized model with the strain fields obtained with the 
full nonlinear model. As previously noted and pointed out in Ref. [H], a Ginzburg-Landau 
type of model can offer an excellent insight even though being defined for infinitesimal dis¬ 
placements. Inspired by these observations, we analyze in detail the linearized model. We 
present the construction of a microstructure which is very similar to both the experimental 
observation and to the numerical solution plotted in Fi g. [T for the full-elasticity model by 
allowing sharp interfaces, that is, when we set K = 0 in (3.1). 

For the reader’s convenience, our construction is split in two parts. In the next subsection, 
we present the geometry of the microstructure. Subsequently, we describe the requirements 
that a map u has to fulfill in order to produce a microstructure which is consistent with 
the experimental observation (see Fig. [^(b)). Theorem]^ which is the core of this paper, 
contains the detailed construction of u. 


3.1 Geometry of the microstructure 

We assume as our reference configuration the set fl := Ri^(0,0), the disk with center in the 
origin and of radius equal to R. We define R := -I- ^ where L is a positive constant. 

In turn, L can be regarded as the characteristic length of the system. To reproduce the 
geometry of Fig. [^we combine a family fcyte-shaped rhomboids that cover (see Fig. 
LEFT). Let t := tan(7r/12) and k G {0} U N. We define the vertices of the rhomboids (here 
and in what follows, when k = 0 we identify Aq = A and similarly for B, C, D, E, F): 


A = 

(xa^Va) 

~ 2 + V 3 ) 

B = 

(xb^Vb) 

_ T(1-V3 1-V3- 

“ 2 V 3 ’ 2 V 3 • 

C = 

ixc,yc) 

_ t(V3 1 I'l 

“ 6 ’ V3 2 I 

D = 

(xD^yo) 

_ rfl 1 1 

“ ^^2 V3’ 2^3 

E = 

{XE,yE) 

_ TI-2-V3 1 ' 

“ 2 V 3 ’ 2V3- 

F = 

[xF^yp) 

_ r/l-|-V3 1 + V3' 

“ 2 V 3 ’ 2 V 3 ' 


Afe 

{xAf^ j yAk ) 

= t^^A 

Bk 

= {xbu.VBu) 

= t'^^B 

Ck 

= ixc^.yc^) 

= 

Dk 

= {xd^iVd^) 

= 

Ek 

= {xE^.yPk) 

= e'^E 

Fk 

= {XFu.yFu) 

= e’^F 


In table we report the complete description of the segments that define the sides of the 
polygons. We can now define the family of open sets whose union gives (up to a set of zero 
Lebesgue measure): 


• for /c = 0: 


IX!A := < 

y{.x,y) G fl : 

y < 9E^{x),x < g^^{y),y> h^^{x),y > h(f'^(a:)|, 


LOb ■ = 

|(x,y) G n : y > g^^{x),y > 3^^(x)|, 


LUc ■ = 

|(x,y) G fl : y > g^^{x),x> y^^(y)|. 


UJd ■ = 

|(x,y) G n : y < g^^{x),y < y^^(x)|. 

xe := -j 

^(a;,y) G : 

y > g§^{x),y < go^{x),y < h^^{x),x < /i^^(y)|. 

Xe := ■{ 

''^{x,y) G n : 

y < g^^{x),y> g^^{x),y < h^^{x),y < /i^^(x)|, 
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Interface 

Equation of the interface 

Bk-A-k 

y = 9k := V^x + xb^ <x < Xa^ 

Ck-^k 

yc,<y< VA, 

CkFk 

y = gk^{x) := -V^x + xc^ < x < xp^ 

DkFk 

y = 9k^\x)-=-^X-^t^^, XD^<X<Xp^ 

EkDk 

y = Qk'^ix) ■=XE,<X<XD, 

EkBk 

y = 9k^ix) := ^x +XE^<x<XB^ 

BkAk+i 

y = ■■=--^x + (^- xb,<x< xa,+. 

Afe+iCfe 

y = := (;^ - < x < xc^ 


y = - <x<xc^ 

DkFk+i 

y = hk^{x) := V^x -b (1 - XD^<x< XB.+i 

DkEk+i 

X = := (I - Vo, < y < 

BkEk+i 

y = hk^{x) := -V3x -b (1 - xb^<x< xb,+i 


Table 1: Sides of the polygons. 


• for /c > 1: 


^Ak ■■= \ 

[{x,y) e : 

y < 9k^{x),x< g^^{y),y> h^^{x),y> h^'^(x)|, 

:= 

|(x,y) e n 

: y > 9k'^{x),y> gk^{x),y < h^\{x),y> hffij. 

^Ck ■= { 

{x,y) e n : 

y > 9k^{x),x > g^^{y),y< h^^^{x),y> hf5L(x)|, 

^Dk ■= { 

{x,y) e fl : 

y < 9k^{x),y < gk^{x),x > h^f^{x),y > /i(?_fi(x)|, 

^Ek ■= { 

(x,y) G n : 

y > 9k^{x),y < g§^{x),y < h^^{x),x < /if^(y)|, 

^Ek ■= { 

{x,y) G fl : 

y < 9 k^{x),y > gk^{x),y < h^^'{x),y< h^^(x)|, 


With the exception of wsjWc and w/j, the sets defined above are rhomboids (see Fig. 
and|^. All these rhomboids are similar to wai in the sense that they can be obtained from 
rescaling and rotation of loa- Indeed, for fc > 1, the sets ujA^T^Ek nnd uiPk are obtained by 
contracting each side of oja,oje and oup respectively, by a factor < 1. Moreover, by a 
^ clockwise rotation of uja we obtain uip and by applying another ^ clockwise rotation we 
obtain ojp. A similar relation holds for ujB^,u}Ck and which are obtained by contracting 
each side of and respectively by a factor of (here with k>2). Then, by 

a ^ clockwise rotation of ujb-i we obtain wci and by applying another ^ clockwise rotation 
we obtain . Finally, by applying a counter-clockwise rotation of ^ and then a contraction 
of each side by t to oja we obtain The measure of the four internal angles of oja are 
respectively ^, f, f, f ■ By similarity, the internal angles of all the rhomboids defined above 
coincide with those of uja- From the equations of the segments reported in Tableit is also 
possible to compute the normals to each side of the rhomboid. By denoting with the 
(two-dimensional) Lebesgue measure, it follows that 

(wAfc) = (ws,) = C? (wb, ) and C'^ (wb^ ) = C'^ (wc,) = (wb, ) Vfc > 0. (3.3) 

Since for A: > 1 wb^ is obtained (up to a rotation) by contracting uJAk-i by a factor t < 1, we 
have C'^iiOBk) = and, since lua^, is obtained by contracting each side of uja^-i 
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by a factor we have £^( 0 ;^*.) = Summarizing: 

L\u:a,) = = LH^% C^cob,) = C^M = = LH^h-\ (3.4) 


If we define 


UJk ■= ^Ak U iOBk U WCfc u U>Dk U UJEk U LOPk 


we have IJ Wfe C fl. Importantly, as the remainder in this inclusion is represented by elements 
of £^-measure equal to zero we have 


= (3.5) 

fc =0 

Notice that, for any fc > 0, the origin does not belong to any of the sets ojk- Furthermore, 
each neighborhood of the origin contains an infinite number of sets Wfc. This fact has deep 
consequence on the property of the microstructure and it is at the origin of the nesting of 
variants close to the origin. 




Figure 3: LEFT: fl and ujq. RIGHT: close-up on uja^, ^Ek k>l. Notice that the 

points Bk-i, Ck-i and Dk-i are not shown in the picture. 


3.2 Construction of the microstructure 

As typical when dealing with microstructure, the definition of the map u is based on the 
introduction of a piecewise-constant tensor field H(a;,?/) such that \7u{x,y) = fi{x,y). It is 
widely known that not every matrix field H(a;, y) is the gradient of a continuous vector field. 
In fact, this condition reveals as a geometrical compatibility constraint on H(a;, y) which is 
known as Hadamard jump condition. Simply, consider a smooth line S having normal v. 
If u is continuous at both sides of S with limits Vu+(a;, y) = H+ and Vu“(a:,y) = H“ at 
each point (x, y) G S', from above and below S respectively, then equating the tangential 
derivatives leads to the Hadamard jump condition 

H+-H-=a(g)n (3.6) 
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Figure 4: ujk, k > 1. 


with a, n e and |n| = 1. The matrices H+ and H in (3.6 1 are said to be rank-one 


connected. Notice that for 2x2 matrices this is equivalent to the condition det(H+—H“) = 0. 
The situation in which we have a map u with Vu piecewise constant and equal to either H+ 
or H- in alternating bands bounded by lines perpendicular to n is relevant in materials 
science. Such maps are called simple laminates m. [a- With no further hypotheses on the 
domain, condition (3.6) is a necessary condition for the existence of a simple laminate u such 
that Vu = H+,H- '. In the case in which we have more interfaces, condition (|3.6|) has to be 


verified at each interface across which we have two crystal phases corresponding to different 
deformation gradients. 

Besides the Hadamard jump condition, which can be regarded as a geometric compatibil¬ 
ity condition, the construction of the map H(a:, y) has to satisfy an energetic condition. In 
fact, this microstructure is zero-energy in the following sense. Recall that £2 and £3 are the 
components of the deviatoric part of E and £1 is the unique component of the spherical part 


the three strains (2.13) then 


of E introduced in (2.11). If we assume that ipL has three distinct minima corresponding to 


/ tpL{£i,£2,£3)dxdy = 0 

Jn 


(3.7) 


if and only if 


{'^u{x,y))^y^ e {Ei,E 2 ,E 3 } for a.e. {x,y) S ft. 


(3.8) 


Therefore, the differential inclusion (3.8) guarantees that the map u is a minimizer for 
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jQfpL{£ 2 ,S 3 )dxdy. Since 'Rsym has to be equal to either Ei, E 2 or E 3 for almost every 
(x, y) S n, what is left to determine is the skew part of H(x, y). 

For the reader’s convenience, we illustrate with an example how the map H can fulfill the 
geometric constraint and the energetic constraint simultaneously. Here we make use of some 
notation that will be further dehned and used extensively in the next section. Consider the 
situation of Fig. [^LEFT, and more precisely the matching between the domains loa and ujb 
. We construct a piecewise constant tensor field H dehned as 


f on UJA 

\ Hb on UJB 


(3.9) 


where the matrices and Hb are to be determined in a way such that {lrlA)sym = E 3 and 
(HB)sym = El and (H^ — Hb) = for some with \v\ = 1. We recall that 

El and E 3 are dehned in ( 2.13| ). Let us write JIa = E 3 and Hb = Ei + W. Here W is an 
unknown skew-symmetric matrix written in the form, 


W = 


0 —w 

w 0 


w G 


(3.10) 


The condition (3.6) applied to and Hb reads (E 3 —Ei — W) = a(g)n or, componentwise. 



= OlTll 

{-^+w)e 

= 01772 


= 02771 

h 

= 02772 


(3.11) 


with nf + n 2 = 1. The system of equations (3.11) has several solutions. Let us choose the 
solution with nf = ^, w = \/3 and vectors a. := (oi, 02 ) = e('\/3, 3), v := ( 711 , 712 ) = (—^, 5 )- 
The vector v is perpendicular to HH, the interface separating uja and ujb- Thus, we obtain. 


Hb = El -|- e 


0 

V3 


-V3 

0 


(3.12) 


The fact that the vector v is exactly orthogonal to the interface separating uja and ujb 
guarantees the existence of a continuous map u such that its gradient coincides with H on 
WA UwbQ 

We now focus on the matching between the domains uja and wbi ■ As the strains associated 
with the domains ujb and wbi are equal, the matching equations to be solved are the same 
than for the matching between uja and ujb, that is, Eqs. (3.11). In this case, however, 
we consider a different solution, namely, nf = w = — ^3 and vectors /3 := ( 01 , 02 ) = 
e(—3, -s/S), M := ini,n 2 ) = (|, ^), where the vector p, is perpendicular to BiA, the interface 
separating uja and wbi- Thus, in this case we obtain. 


Hr, =Ei 


0 

-^3 0 


(3.13) 


Notice that the symmetric part of Hr^ (as well as the symmetric part of Hr as defined in 
Eq. (3.12)) is precisely Ei. The fact that y, is orthogonal to the interface separating uja and 


wrj is enough to guarantee that there exists a function u such that Vu = Ha on uja and 


^We note that the geometric constraint has been implicitly taken into account in the construction of the 
star-shaped geometry of the previous section. 
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Vu — on uiBi- Of course, if u = u over uja one can define a continuous function defined 
over the larger set uja U wb U wbi ■ 

Even though the example above is restricted to three subsets of O, it suggests that this 
technique requires a careful matching of matrices and depends heavily on the geometry of 
the problem. Whether it is possible to further iterate this construction over a larger domain 
is not a trivial question. Here one of the mathematical difficulties lies in the fact that this 
microstructure becomes finer and finer close to the center of the tripole-star pattern. Indeed, 
as fc —> 00 , the size of the sets ujk is reduced while the number of interfaces increases. Since 
one must be able to prove that the Hadamard jump condition is satisfied at each interface, it is 
not a priori clear if the construction of a continuous map u can be extended asymptotically 
for each hierarchy of laminates. In other words, the point (0, 0) may act as a topological 
singularity. Below we show that this construction can actually be extended to the whole 
n. Additionally, in accordance with experimental observations, our construction involves a 
mixture of variants occurring with the same phase fraction. By defining 

WE. := |(a;,y) G ; (Vu(a;,j/))^^^ = E,|, i = 1,2,3, (3.14) 


the sets of points {x, y) in B such that (Vu)sym = E^, we have that wei, wej; and wes have 
the same area. 

In the remaining part of this section we construct the function u from scratch instead 
of focusing on the abstract properties for which such a construction can be achieved. This 
approach has, of course, an intrinsic limitation in the sense that it is not clear if it can be 
extended to treat more general cases of non-periodic microstructures. The construction of u 
will serve as the proof of the following theorem, which is the central result of this paper. 


Theorem 1. Let L and e be any two p ositive real numbers. L et VL = 
-I- let Ei be defined as in (2.13) and we; as defined in (3.14). 
map u : B —>■ such that 


Br{ 0,0) with R = 
Then there exists a 


1 ) \i is continuous in Q, 

2) {Wu{x,y))^^^ G {Ei,E 2 ,E 3 } for a.e. (x,y) inn, 

3) C^(ujE,) = ^R^ * = 1,2,3, 

4) u lies in (an equivalence class in) IT^’^(B,K^), 1 < p < oo. More precisely, u is not 
Lipschitz continuous. 


Remark 1. Proof of Theorem [7] lies in a tedious but straightforward series of algebraic 
computations. To keep the presentation engaging, we give here a sketch of the argument 
required. To construct the function u, we must ensure that the geometrical compatibility 
condition holds at each interface. Indeed, consider the two rank-one connected matrices Ha 
andtlB, introduced in the example above. Since compatibility holds at the interface BA, upon 
integration of Ha andHB we can define, up to constants of integration, the two functions 


evA{x,y) 



“I “ + '^^.1 

— X^ -f I -f CA,2 


evB(a;,y) 


X - V3y cb,i 
V3x -y-\- cb ,2 


over u!A andujB respectively. By requiring continuity at the interface we determine cb,i and 
cb, 2 , thus leaving ca^ and ca ,2 unknown. Similarly, matrices He, Hr, He andHp can be 
defined by requiring their symmetric part be equal to either Ei,E 2 or E 3 and by requiring 
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them to be pairwise rank-one connected thus fixing their skew parts. Upon integration of 
these matrices one defines vectors vc> v_d, ve andvp overujc, <^d, <^e andujp respectively. 
Notice that the constants of integration will again he determined by matching the vectors at 
the interfaces (except ca,i and ca ,2 which are undetermined). We have then constructed a 
piece-wise affine vector field which is continuous over loq and which we denote with Vq. An 
analogous procedure over ujk yields a piecewise affine vector field which we denote with Vk. 
As for Vo, each function with k > 1 is defined up to constants of integration. Finally, the 
introduction of the globally continuous map u follows after the computation of the remaining 
constants of integration (up to two which remain, correctly, unknown) by matching the maps 
"Vk and Vfc_|_i for every k > 0. This is in fact possible because the geometric compatibility 
condition holds also for each interfaces between uJk and cok+i ■ 


Proof of Theorem [TJ In what follows we make an extensive use of the notation introduced 
in Paragraph |3.1[ Let us define 


where v : IJ^g 


u(a;, y) = ev{x, y) + e\o{x, y), 
is defined as v = v^ on Wfc and v^ : ujk 


(3.15) 


is defined as 


Vk{x,y) := < 


VAfc on WAfc 
vsfc on ws, 
vc, on wc, 
VDk on LUDk 
ve^ on uje ,, 
VEfc on WEfc 


fc > 0 


(3.16) 


with 


^Adx,y) ■= 


^Bk{x,y) ■■= 


^Ckix,y) := 


'^D„{x,y) := 


_| + (2fc-i)y3y 

-V^x(-+2fc) +1 


Lt 


2k 


j.2k 


X + {—Vs + 2kVS)y + 


4.2k 


x{VS - 2kVS) - 2 / + ^(2 


VS)L 


X 

"2 + 


VSy(^ 


2k-- 

2 


_ 2kV3x + I 


+ 2kVSy - 


sVs 






(^ _ 
A 2 


2 kVS )x + - + 


(A^) 
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^EAx,y) ■■= \ ^ 1 AJ 

[ (- - 2 fc ) + I + + _)l 

{ X + 2V3ky — 

/o 

—2kV3x — y + 

and Vo : IJ^o defined as Vq = Vofc on Wfe and Vofc : Wfc —)■ is defined as 


' fc 

L(V^ - 2) ^ 

i=i 

. j=i 

By construction, the function u is defined at almost every point in Q. First of all, we show 
that there exists a continuous extension u of u defined at each point {x,y) of fl. Notice that, 
on each subset ujAk > ^Cu ■> ^Dk > “^-Efc the map u is affine and therefore continuous. 

We now check that u has indeed continuous limit across each interface. In what follows we 
compute the limit u+,u“ at each interface. 

To begin, consider all the interfaces inside Wfe. Let fc > 0: 

1. Interface B^Ak 



for fc = 0 : Voo(a;,y) := 


for fc > 1 : ^ok{x,y) ■= 


u {x,y):= lim \i{z^,Zy)= lim e{vB^+^rok){zx,Zy) = 

{z^,Zy)^{x,y)- {Za:,Zy)^{x,y) 

= lim t{^rAk+^ok){Zx,Zy)= lim u(a;,y) =: u+(a;,y) = 

{z^,Zy)^{x,y) ( 2 ^, 2 y)-^(a:,y) + 

r —2x + (!)kx + L{2k—\)t^^ 
= e^ok + e I -2y/ikx + 


2. Interface EkBk 


lim 

{z^,Zy)^{x,y) 


+'Vok){Zx,Zy) = 


lim 

{z^,Zy)^ioi:,v) 


= eVofc + e■ 


eCvBfc +'Vok)izx,Zy) = 

' 2fcx + Lt2fc^+2.t2fc^(V3_2) 


2^3 


^x — 2k\/ix + 


73 


3. Interface EkDk 
lim 


, , , A^Du+^ok){Zx,Zy)= lim AEu+^ok){Zx,Zy) = 

{Za:,Zy)^{x,y) {Zj:,Zy)^{x,y) 


= evok + e ■ 


LAk-A{^ - \)L 

^x - 2kV5x + 
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4. Interface D^Fk 

lim e(v£,, + ^ok){zx, Zy) = lim e{vF^ + y^ok)izx, Zy) = 

(z^,Zy)^{x,y) {z^,Zy)^{x,y) 

j X — 2kx — 

= ev„, + e| 


5. Interface CkFk 

lim t{^rck + ^ok){zx,Zy) = lim e{^^F^ + ^ok){zx, Zy) = 

{z^,Zy)^{x,y) {z^,Zy)->-{x,y) 


= eVok + e ■ 


X — 6kx — ^2 -1" 

—2ky/ix + VSx + 


6. Interface CkAk 

lim e{^^c„ + ^ok){zx, Zy) = lim e(vA, + ^fok){zx, Zy) = 

(■2x,2y)-^(ai,y) (2x,2j,)-)-(a:,y) 

_ , / + 2k\/iy - 


Let us now consider the interfaces between w/c and Wfe+i, for any k > 0: 


1. Interface 

u"(x,y):= lim u(z^,Zy)= lim e(vB,+i + Vofc+i)(z^, = 

{Za:,Zy)^{x,y)- {Za:,Zy)^{x,y) 

= lim d-VA^+-Vok){zx,Zy) = lim u(a;,y) =: u+(a;,y) = 

{z^,Zy)^{x,y) (2x,2y)-^(a:,y) + 

-2kx + 2kLf'^{^ - 1 ) + e’^L{\ - 
- 2'/ikx + 


= evofe + e ■ 


'3 2^/3' 


2. Interface B^Ek+i 

lim e(vB;,+i + ^ok+i){zx, Zy) = lim e(vB^ + Vok){zx, Zy) = 

{z^,Zy)^{x,y) {Za;,Zy)^{x,y) 

j —2x — 6kx + Lt'^^{—k — 4:k + 2kVS) 
= + M -2fcx^/3 + L^2fe(I + ^) 


3. Interface D^Ef^Fi 

lim + Vofc+i)(z^,Zj;) = lim +^ok){zx,Zy) = 

(2x,2y)—>-(ai,y) {Za:,Zy)^{x,y) 


= + e ■ 


v/3(i+2fc)y + Lt2fc(i-^) 

f+Lt2fc(2fc_A:V3+|V3) 


4. Interface DkFk+i 

lim e{vD,,+,+Vok+i){zx,Zy)= lim e(vF^ + Vok)(zx, Zy) = 

{Za:,Zy)^{x,y) {Za;,Zy)^{x,y) 


= + e « 


x + 6kx — + 2k\/iLt^*^ — AkLt'^'^ 

x{-V3 - 2fc^/3) + Lt 2 fc(_i + 
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5. Interface Fk+iCk 


lim 


lim +Vok+l)iZa:,Zy) = 

{z^,Zy)^{x,y) (z^,Zy)^(x,y) 


= eVok + e ■ 


e(v_Ffe + Vofc)(zx,Zy) = 

x{l + 2k)+Lt^>^{2k-1) 
x{-2kV^ - ;^) + 


6. Interface Ak+iCk 


lim 


>{x,y) 


e(vc',+i +'^ok+i){zx,Zy) = 


lim 


e(vAfc +'Vok){Zx,Zy) = 


= £^ok 


(z^,Zy)->-(x,y) 

r + 1 ^ 2 fc-y 3 fc) 


-^x- 2kV3x + + !)■ 


2 V 3 


The above calculations show it is possible to extend u to a function still denoted by u by 
continuity across each interface. Notice that the function u is continuous also at the vertices 
Ak, Bk, Ck, Dk, Ek and Fk. The new map u is now defined at each point of n/(0,0) and, 
for each 5 > 0, u is continuous on the set 0). Therefore, in order to extend u to a 

continuous function on the whole il, we have to clarify the behavior of u close to the point 
(0,0). Precisely, we prove that 


lim u{zx,Zy)=eL{V3-2,l)-^—;z. 

{Z^,Zy)^{0,0) 


(3.17) 


Notice that it is enough to show that the supremum and the infimum of the components of u 
on LUk generate sequences {M^} = {Mi k,M 2 ^k} and {mj.} = {mi^k,'m 2 ,k} both converging 
to eL(v^- 2,1)^ as fc —> 00 . This, combined with a standard contradiction argument 
guarantees that (3.17) holds for every sequence of vectors converging to (0,0). Let us denote 


with Ui(x^y), with i = 1,2, the components of u(a;,j/). We compute the supremum of Ui by 
splitting the computation over the subdomains composing ujk yielding 


Mi^k ■= sup Ui{x,y) =ma.xU sup Ui{x,y)y( sup Ui{x,y)\, 

{x,y)^uik ^ ^{x,y)^u]A^. ^ ^{x,y)euiB^. 

{ sup Ui{x,y)y( sup Ui(x,y)],f sup u^(x,y)],f sup Mi(a;,?/)j|. (3.18) 


^(x,y)euiD^ 


^{x,y)eLLiF^ 


Since u is affine over each of the subdomains in (3.18), supremum and infimum are attained 


on the boundary (and, in particular, on the vertices). Therefore, focusing on the case of wa*. 
we can write 

sup u^{x,y) = max{zq(xAfc,2/Aj, Ui{xB^,yBk), u^{xCk,yCk), Mi(a:Afc+i,yA^+i)}-(3.19) 

{x,y)^u]Au 


A straightforward computation yields 

^ I Lt- [^ + 1 - ^] ELY.U 


(3.20) 


u{xBk,yBk) = £■ 


[fc(- 1) + i - ^] + (v^ - 2) L Y!;=i 
[ k{V3 - 1) + 1 + ^] + i E ■=! 


-2 


(3.21) 
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(3.22) 


u(a:Cfc,2/cJ = 

u(xAfc+i,?/Afc+i) =e| 

Recalling that 


f Lt^/^[k(2-V3) + ^-l] +(V3-2)LZj^it^^-^ 

i + (fc + l)(2 + y3)] + (V3-2)LZj=it^^-^ 


+ 1 - (fc + 1)] + LE-i 


(3.23) 




i=i 


-t2 


it is now immediate to observe that u{xA^,yAk)^ ^{^CkTVCk) 

^{xAk^nDAk+i) converge to eL(-\/3 — 2 , 1)1372 as k ^ 00 . Repeating this argument 
for the remaining regions lOBk-, ■ ■ ■ ^^Fk yields 


Mfe = (Mi,fc,M 2 ,fe) eL(y3- 2 , 1 )--- 5 - as k ^ 00 . 

1 — 

The, we can define 

rrii^k ■= inf Ui{x, y) = mini ( inf Ui{x,y)),( inf Ui{x,y)), 

(x,y)GuJk ^ \(x,y)GuJA,_ / \(x,y)Giyn,^ ' 


(3.24) 


( inf Ui{x,y)),( inf Ui{x,y)),( 
\{x,y)GujCk ' \(a;.iy)ecJr). / V 


inf Ui 
{x,y)^y)Dk 


inf Ui 

^(x,y)eyJE^ 


(x,y)eujBk 

i(a;,y)),( inf M;(a;,?/))|. (3.25) 

Recalling that the infimum of u is attained on the vertices, the same computation of 
(3.20 3.23) is enough to prove that the sequence {inE^^ (ui), inE^^ (uj)} converges to 
eL{\/3 — 2,1) as fc —>■ 00 . Repeating the same argument for the remaining infima in 
(3.25) shows that 


mfc = (TOi.fc,m 2 ,fc)-» eL(v^- 2 , 1 ) ^ as fc 00 . 

I — 

Finally, by combining (3.24|) and (|3.26[) we have 


u{x,y) -A eL{V3- 2,1)- 


1 


as (x,y)-^(0,0) 


(3.26) 


(3.27) 


l-f2 

as claimed. Therefore, by a further extension of u we can define the continuous function 

f u(a;,j/) on fl/( 0 , 0 ) 

u(x,y):=<| eL(v^-2,l)Y^ _ for (x, y) = (0,0) (3.28) 

[ extended by continuity on dfl. 

The fact that u^(xA,yA) = u_b(xa, 2 /a) = uc{xA,yA), uc{xF,yF) = u_f(xf,?/f) = 
UD{xF,yF), and nF){xE,yE) = ^E{xE,yE) = ^BixEjys) guarantees that the function de¬ 
fined above is indeed continuous up to the boundary of ft. 

We now prove point 2). Since u(x,?/) is affine over each of the sets ujAk, ^Bk^ ^Ck^ '^Dk^ 
ujEk and uJFk, h is straightforward to compute its gradient: 


Vu(x, 2 /) = H(x, 2 /) := < 


llAk{x,y) 

on 

UJA 

Hsfc(x,y) 

on 

ojb 

^Ckix.y) 

on 

UJC 

^Dkix,y) 

on 

lud 

^Ek {x,y) 

on 

uje 

. '^Fk{x,y) 

on 

OJFf 


(3.29) 
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where 


-V3{2k+l) 


y3(2fc- i) 

1 

2J 2 


2' ) =E3+Wfc 


Hr, =e 


{-V^ + 2kV^) 


(V3-2fcy3) -1 


= El + Wfe + W 


(3.30) 


(3.31) 


He. = 


^-2ky/3 


V3(2fc- i) 


= E, + Wfc + W 


(3.32) 


He. =e 


2ky/l-^ 


{^-2kVS) 


= Ea + Wfc + W 


(3.33) 


and 


He. - e 


(i-2A:)V3 


{}+2k)V3 \ 

I J 


— E 2 + Wf. 


Hf. =e 


2A:V3 


-2kV^ -1 


= El + Wfe 


w I 0 2^73 ^ / 0 -v^ 

'^* = '1 -2.v^ 0 )■ ^ = ‘[ V3 0 


(3.34) 


(3.35) 


(3.36) 


Therefore we have (Vu(a;, y))^^^ G {Ei,E 2 ,E 3 } for almost every {x,y) in Q. 

Proof of point 3) is also straightforward. By combining (3.29) with ( 3.30|3.35 1 we can 
write 


C'^(^{x,y) G n : (Vu(a;,y))^^^ = Ei) = +£^{1^3,)^ 

k^O 

00 

£^(^{x,y) : (Vu(a;,y) G = E2) = '^(£^{u}e^ +£^{u}cJ), 

k^O 

00 

£^(^{x,y) : (Vu(a;,y) G = E3) = +^^(we.))- 


k=0 


Then, recalling that by (3.3) we have C‘^{ujAk) = ^^(“^F.) = ^^{^Fk) and C’^iujE^) = 
£^(a;efc) = £^{^Dk) for every fc > 0, we conclude that 

Y2(^£'^{ujFk+£'^{uJBk)j =Y^(£'^{uje^+£'^{ujcJ^ = ^(^/:^(a;A.+/i^(a;e.)) = —■ 


fc=0 


fc=0 


fc=0 


We are left with the proof of point 4), that is, u G with 1 < p < oo. We only 

have to show that Vu G with 1 < p < oo. Componentwise, we have to show 

that Vuij G L^(f2) with = 1,2. For {x,y) G ujk we have 

\Vu,j{x,y)\P < 2P~'^{\\'Usym,ij{x,y)\P + NUskew,zjix,y)\P). 
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Notice that by (3.30)-(3.35) ^Uaym,ij G L°°(Vt) for i,j = 1,2. Therefore, to obtain the claim 
we only have to show that Vuskew,ij is in with 1 < p < oo. Of course, it is enough to 

consider the case i ^ j. By Beppo Levi theorem on monotone convergence we have 


/ \{^^skew')ij\^dx(ly — / 'y ^ \{^ Uskew')ij\^\iVk(.^: y') — 

Jn JQ 

N „ ^ r 

lim V / |(Vu skew)ij\^XujAXyy)dxdy = lim ^ / \{Vuakew)^ 3 \^XuJki^^y) dxdy + const., 


A ;=0 ‘ 


k^l * 


where the characteristic function Xujk equal to 1 on ujk and 0 in Q\ujk- For A: > 1, in view 
of (3.30)-(3.35), we can write 

clk^Xi^ki.x,y) < \{'^uskew)t3?xu,ki.^,y) < clk'^Xujki.x.y)- (3-37) 

Here ci and C 2 are two (positive) real numbers that do not depend on k. Therefore we have 

w » N „ 

lim V / \(yuskew)ij\^Xujt:ix,y)dxdy < lim W / dxdy. (3.38) 

AT-foo Af-i-oo ^ 


k=l 


k=l ' 


To conclude it is enough to estimate the right hand side in (3.38). Recalling (3.4) we can 
write 


fc=i' 


lim f k'^’XbJk y) dxdy = ^im 3 (^j k^ dxdy + f kP dxdy^ = 

N 

lim 3 (3.39) 


k=l 


□ 


Since t = tan(7r/12) < 1 this series converges for 1 < p < 00 . 

Remark 2. Notice that from point 4) of Theorem^ it automatically follows by the Sobolev 
imbedding Q Thm. 6.1-3] that u S C'®’'’(H,R^) H (7(11, R^), where C'°’’'(H, R^) is the space 
of Holder continuous functions with 0 < r < 1. However, notice that our proof of point 4) 
implicitly makes use of the fact that u does not jump across the interfaces, which follows 
from 1). Moreover, our proof of 1) has the advantage of actually determining the value of 
the map u in the origin. 

Interestingly, from Point 4) of Theor em\]\ we conclude that u is not Lipschitz continuous 
in XI. To verify this, observe that by Eq. [3.31) we can not find an uniform bound for X/Uskew 
in Xt. This is a relevant difference with respect to the theory of simple laminates m, m, 
m in which solutions of differential inclusions are usually sought in the space of piecewise 
affine amps. However, given any positive 5, u is Lipschitz continuous over the smaller set 
Xl\Bs{0,0). 


4 Discussion 

Theorem [l] guarantees the existence of the function u satisfying the properties stated in 
points 1-4) of the theorem, but it does not guarantee its uniqueness. The absence of addi¬ 
tional constraints in the model, such as boundary conditions or applied forces, enables the 
construction of a different map u still satisfying all the points of the theorem possibly with 
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a different geometry. For instance, the definition ( |3.28 ) can be modified by adding a rigid 
deformation (that is, a map whose gradient is a constant skew-symmetric matrix) and the 
new function u still fulfills points 1-4) of Theorem[^ Indeed, let z = (zi, Z 2 , ^a) be any vector 
in and let 


u^(x,y) = u{x,y) + 


0 

-2:1 


^1 

0 


^2 

Z3 


(4.1) 


where u is defined in (3.28). Then, Uz(x, y) clearly fulfills points 1-4) of Theorem]^ Letting 
0 = z = (3^,0,0), in Fig. ^we plot the graph of the function VLz{x,y) and in Fig. ^ we 
plot the level curves of u^. The choice of z is purely for clear representation purposes. From 
these figures it becomes clear that the field \iz{x, y) is locally piecewise affine (except at each 
neighborhood of the origin), and that there is no singularity at the center of the tripole-star 
pattern. In particular, it has to be noticed that the displacement \iz{x, y) is continuous 
across the interfaces. This property is a direct consequence of the geometrical compatibility 
condition which has been imposed in the construction of u{x,y). While the construction 
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Figure 5: Plot of u(x,y) = {uzi{x,y),Uz 2 {x,y)). LEFT: Uzi{x,y). RIGHT: Uz 2 {x,y) (here 
for simplicity L=e=l). The web version of this article contains the above plot hgures in color. 




Figure 6: Level curves of U 2 (x, y) = (■U 2 i(x, y), Ziz 2 (a:, 2 /)). LEFT: U 2 i(x, 1 /). RIGHT: ^^ 2 (a:, y) 
(here for simplicity L=e=l). 

of Theorem works for any pair L > 0, e > 0, it is clear that some limitations have to 
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be imposed on L, e in order that the assumptions of linearized elasticity are fulfilled. Since 
u is continuous in it is straightforward to derive conditions on L, e that guarantee that 
|u| = o(l). In fact, by denoting with ui(x,y) and U 2 {x,y) the components of \i(x,y), we can 
estimate _ 


max_ |u(a;, j/)| < 

(x,y)Gfl 


max_|Mi(a;,y)| 

L(a;,y)GO 


2 

+ 


max_ 

{x,y)en 


\u2{x^y)\ 


2 


where 

max_ j/)| = max< yjia.x_Ui{x, y), — min_Ui(a:, y) k i = l,2. 

{x,y)£Q ^{x^y)GQ, (x,y)GQ, ^ 


As in the proof of Point 1) of Theorem]^ we are left with splitting the computation of the 
maximum and the minimum of the components of u over into the evaluation of ui{x,y) 
and U 2 {x,y) on the vertices of the rhomboids: 


max_ u,(x,y) = max<^ u,{xAk,yAj,Ui{xB„,yBk), ,yc^),Ui{xD^,yDj, 

(x,y)en '■ 

Ui{xEk,yEk),Ui{xFk,yFk),u^{0,0), A:e{0}UN|, i = 1,2,(4.2) 


min_ Uj {x,y)= min<^ u* {xa, ,yAj, Ui (xs,, ys J, (xc,, yc J, m* {xd^ ^yoj, 

(x,y)en '■ 

UiixE^,yEi,),Ui{xFk,yFk),u^{0,0), A:e{0}UN|, J = 1,2.(4.3) 


Now, by simply plugging the equations of the vertices (3.2) into (3.28) we immediately obtain 
that 


max_Mj(x, y) oc eL. (4-4) 

(x,v)Gfl 


Therefore, if eL ^ 1 we have that |u(x,y)| = o(l) for every (x,y) in Q. 

The discussion for the gradient of the displacement is more delicate. The fact that Wugkew 
is unbounded in 12 (see Remark has deep consequences on the validity range of our model 
in the scenario of linearized elasticity (see 0 Chapter 11). Correctly, such a model holds 
under the assumption that |Vu| is small. This poses a requirement on the symmetric and 
skew symmetric parts of Vu that have to be simultaneously small. While by Eq . (2.14) 
|Vusym(a^,y)| = v^e is uniformly bounded at almost every (x,y) in 12, Eq. (3.37) tells us 
that |Vusfeeu,| grows linearly in k over ujk- Therefore for k large enough, the assumptions 
of geometrically linearized elasticity fail. Whether replacing the linearized theory with the 
full non-linear elasticity model would bring more insight onto the analysis close to the origin 
is not clear. Although the full analysis of nested microstructure of the type of Fig[^ in the 
more general model of finite elasticity is to date an open question, here we try to address 
this point with an example. Consider for instance the point 


1 / 5 5 \ 

Cfe := (x^k^y^k) = +a;B,+i,yBfc +2 /b,+i) -3,3- 

which belongs to lob^ and is equally distant from and Bk+i- As fc —>■ oo the sequence 
{?fc} converges to (0,0) at a rate of . In the case where L « 0.07ym (as can be deduced 
from Fig. [^(a)) we have that, already for /c = 2 

1 ^ 2 ! 10“^nm 
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which is a distance smaiier than the crystai iattice parameters. Thus, in order to capture 
aii the detaiis of the behavior of the system ciose to the origin, even the use of a fuii non- 
iinear eiastic modei wouid not be justified. This is where atomistic modeis wouid be more 
reiiabie. Nevertheiess, the domains obtained for fc = 0,1 and corresponding to the regions 
Wq = loa U wb U ujc Uwd U uje U cop and ujb-i^^Ci and remain in the range for which 
continuum approximations are acceptabie. Notice that due to the fast convergence in fc, 
these regions cover the great majority of our reference domain fl. Indeed, a straightforward 
computation shows that the area of this region corresponds to 




0.999, 


thus confirming that the region corresponding to fc > 2 is confined in a tiny neighborhood of 
the origin. 

We want to remark that the domain-length L = 0.07/rm reported above represents only 
a possible example of an observed microstructure. Due to self-similarity, other tripole-stars 
with larger L may well be observed for which higher fc-order generations remain in the range 
where continuum models are valid. 

Finally, we discuss the effect of neglecting the interfacial energy. Computations based on a 
model for the triangle-to-centered rectangle transformation with sharp interfaces (A. Ruland, 
forthcoming) show that minimizers of the total energy may develop singular gradients also in 
the presence of interfacial energy terms. Formation of the self-similar microstructure is not 
prevented by an energy penalisation of the interfaces because the line energy contribution 
associated with the nesting is negligible. This energy argument thus validates our findings 
based on a model with no interface energy contribution. 


In summary, we have analyzed a nested martensitic microstructure observed in lead or¬ 
thovanadate and Mg-Cd alloys. In the scenario of linearized elasticity, we have modeled the 
microstructure as a three-phase martensitic mixture by showing that geometric compatibility 
across the interfaces holds for each interface between the variants making up the hierarchical 
microstructure. Continuity of the displacement (Theorempoint 1)) rules out irreversible 
phenomena such us cracks and cavitations in the model analyzed in this article. This work 
also offers an analytical counterpart to the numerical work of Ref. m- The approach 
adopted in Ref. |12j is strain-based, in that the minimization problem is written in terms of 
strain components (rather than on the displacement map u) and geometric compatibility is 
imposed through a partial differential equation relating the strain components. Thus, the 
properties of u are not investigated. In Theorem 1, we have proved that geometric compati¬ 
bility is not only achieved globally in the domain but also close to the origin. We have also 
shown that the map u, which realizes the microstructure, is smooth. 

One of the main contributions of this paper is the precise estimation for the growth rate 
of the gradient of the displacement and on the computation of the geometry of the domains 
close to the singularity. The unboundedness of the gradient of the displacement is an intrinsic 
aspect of the analyzed system. While no continuum model, either in non-linear or geometri¬ 
cally linearized elasticity, is appropriate to exactly describe the singularity, we provide precise 
quantitative information on the scaling law of the self-similar structure and estimate the size 
of the region where a linearized model becomes inconsistent. Remarkably, the geometry of 
the microstructure we compute corresponds to the one observed experimentally and to the 
ones computed numerically with a good level of approximation (see Fig. [^. 
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Whether the results of this paper can be extended to treat more general situations be¬ 
longing to the large set of non-periodic microstructures and, possibly fractals, is an intriguing 
question. This will be considered in a forthcoming paper. However, we have addressed one 
of the most challenging mathematical aspects associated with this question. Nesting of mi¬ 
crostructure close to a singleton may turn into an increasing gradient of the displacement. 
Although in the current situation this results in loss of regularity of u which fails to be Lip- 
schitz in fl), it is not clear if the approach described in this paper may be applied to more 
general situations. Thus, a more general theory of nested martensitic microstructure requires 
treatment of geometric compatibility in the presence of potential point defects. Domains with 
singularities may be seen as special cases of non-simply-connected sets and therefore would 
require the types of techniques developed in Ref. m ( see also Chapter 2]). 
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